Intermittency in the isotropic component of helical and non-helical turbulent flows 
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We analyze the isotropic component of turbulent flows spanning a broad range of Reynolds num- 
bers. The aim is to identify scaling laws and their Reynolds number dependence in flows under 
different mechanical forcings. To this end, we applied an SO(3) decomposition to data stemming 
from direct numerical simulations with spatial resolutions ranging from 64 3 to 1024 3 grid points, 
and studied the scaling of high order moments of the velocity field. The study was carried out for 
two different flows obtained forcing the system with a Taylor-Green vortex or the Arn'old-Beltrami- 
Childress flow. Our results indicate that helicity has no significant impact on the scaling exponents. 
Intermittency effects increase with the Reynolds number in the range of parameters studied, and 
are larger than what can be expected from several models in the literature. Finally, we confirm 
previous results showing the use of the SO (3) decomposition improves the convergence of scaling 
exponents. 

PACS numbers: 47.27.ek; 47.27.Ak; 47.27.Jv; 47.27.Gs 



I. INTRODUCTION 

Turbulence is a recurring phenomenon in nature; we 
can find turbulent dynamics in atmospheric, geophysical 
and astrophysical flows. The dynamics of these flows is 
often described by identifying three characteristic ranges 
of scales the injection range, whose properties depend 
on the forcing; the inertial range, which is assumed to 
have universal properties in the limit of infinite Reynolds 
number; and the viscous range, where dissipation takes 
place. In three dimensional isotropic and homogeneous 
turbulence, the scale separation between the forcing and 
the viscous range (i.e., the width of the inertial range) in- 
creases as a known power of the Reynolds number. Nowa- 
days, computing power is scarcely sufficient to study 
in direct numerical simulations (DNS) flows with these 
three ranges well resolved. Even in the few simulations 
where an incipient scale separation is achieved (USUI, an 
exploration of the parameter space to build confidence on 
such assumptions as universality of inertial range prop- 
erties is currently out of reach. 

In a recent study [|| numerical simulations up to spa- 
tial resolutions of 1024 3 grid points were performed using 
different forcing functions, including coherent and delta- 
correlated in time forcing, as well as using mechanisms 
that injected only energy or both energy and helicity into 
the flow. From numerical simulations [f| 0, H, 0] it is 
known that helical and non-hclical isotropic and homo- 
geneous turbulence follows a Kolmogorov scaling in the 
inertial range, albeit intermittency corrections. However, 
in Q a departure in the scaling exponents of higher or- 
der moments between helical and non-hclical simulations 
was found. It was unclear whether these departures were 
associated to a dependence with the helicity content of 
the flow, or with anisotropics generated by the different 
forcing functions acting on the large scale. Other stud- 
ies reported a dependence of the intermittency of the 



energy flux [3] or of the recovery of isotropy on helicity- 
dependent statistical quantities (To| on the amplitude of 
the helicity flux. 

In this work, we use the SO(3) decomposition [Til. Il2| 
(see also [l3| for a review) to separate the isotropic and 
anisotropic components of a turbulent flow. To deal 
with the size of our datasets, the decomposition is imple- 
mented numerically as described in [f~4| . Previous studies 
focused on the scaling of the anisotropic sectors (see e.g., 
HUGHES]) from data stemming from numerical simula- 
tions or experiments, as it has been conjectured that both 
the isotropic and anisotropic sectors of turbulent flows 
may follow universal scaling laws in the inertial range 
(see e.g., [TU, [H|), an d as f° r the third order structure 
function it has been rigorously shown that the "four-fifth 
law" holds in the isotropic sector even in the presence of 
anisotropies [ll|, Qj|. These studies also showed (albeit in 
simulations or experiments at moderate Reynolds num- 
bers) that the use of the SO (3) decomposition improves 
the convergence of high order moments of the velocity 
field [HI, QjJ ■ F° r the determination of just the third 
order scaling exponent, studies of the isotropic compo- 
nent has been done up to spatial resolutions of 512 3 grid 
points 1 1 -41 ] . However, less attention has been paid to com- 



parisons of the scaling exponents in the isotropic sector 
for orders other than second or third, and to compar- 
isons of these exponents for different forcing functions in 
simulations at large spatial resolution. Considering the 
differences reported in helical and non-hclical turbulence, 
it seems reasonable to consider such flows for a compari- 
son. Here, we use velocity fields stemming from six DNS 
with spatial resolutions ranging from 64 3 to I024 3 grid 
points, using either helical or non-helical forcing mecha- 
nisms. Velocity increments and scaling exponents in the 
inertial range are computed for all runs, up to eight order 
for the simulations with the larger spatial resolution. 

The analysis indicates that helicity has no measurable 
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effect on the scaling exponents of the velocity field, with 
the differences observed in the exponents stemming from 
the different runs being within error bars, or ascribablc 
to Reynolds number dependence. Moreover, we find that 
intermittency (as measured from the departure of the 
linear dependence of the scaling exponents with the or- 
der, associated to the development of strong events in 
the velocity field) increase with the Reynolds number in 
the range of parameters studied, indicating convergence 
of high order moments to their infinite Reynolds num- 
ber limit is slow and may have not been achieved even 
in the largest numerical simulations up to date. Finally, 
we confirm previous results showing that the use of the 
SO (3) decomposition improves scaling laws and the de- 
termination of scaling exponents. 

The structure of the paper is as follows. Section |TT] 
briefly describes the numerical simulations. Section Ifffl 
introduces the SO (3) decomposition and describes the 
method we use to apply the decomposition to our data, 
based on the implementation proposed in [f4j . Section 
IIVI discusses the 4/5 law and the energy scaling in the 
isotropic sector, and Sect. [V] presents the scaling laws 
obtained in the inertial range for the different moments 
of the velocity field. Also, a comparison between helical 
and no-helical flows is done, as well as a comparison with 
models of intermittency. Finally, Sect. I VII presents the 
conclusions. 



II. NUMERICAL SIMULATIONS 

The datascts we use for the analysis stem from DNS 
that solve the momentum equation for an incompress- 
ible fluid with constant mass density. The Navier-Stokcs 
equations under these conditions read 



+ v . V v = fl+f- 

at p 



zA7 2 v, 



(1) 

V-v = 0, (2) 

where v is the velocity field, p is the pressure, f is an 
external force that drives the turbulence, v is the kine- 
matic viscosity, and p is the mass density of the fluid (set 
to unity here in dimensionlcss units). The mode with 
the largest wave vector in the Fourier transform of f is 
defined as the forcing wavenumber kf, with the forcing 
scale given by Lf = 2n/kf. 

Equations {J) and ([2]) are solved using a parallel 
pseudo-spectral code in a three-dimensional box of size 
2tt with periodic boundary conditions [2(| Hl[ . We use 
three different spatial resolutions: 64 3 , 256 3 and I024 3 
grid points. The equations are evolved in time using 
a second order Runge-Kutta method, and the code uses 
the 2/3 rule for dealiasing. Reynolds numbers quoted are 
based on the integral scale and defined as Re = UL/v, 
where U = {v 2 ) 1 / 2 is the r.m.s. velocity and the integral 
scale L is defined as 

/ E{k)k~ 1 dk 



TABLE I: Parameters used in the simulations. N is the linear 
resolution, f is the forcing (either TG or ABC), kf is the 
forcing wave number, v is the kinematic viscosity, and Re is 
the Reynolds number. 



Run 


N 


f 


kf 


V 


Re 


Tl 


64 


TG 


2 


5 x W~ * 


40 


T2 


256 


TG 


2 


2 x 10~ :i 


675 


T3 


1024 


TG 


2 


3 x 10~ 4 


3950 


Al 


64 


ABC 


3 


4 x 10~ 2 


70 


A2 


256 


ABC 


3 


2 x 10~ 3 


820 


A3 


1024 


ABC 


3 


2.5 x 10 -4 


6200 



L = 2tt; 



jE{k)dk 



(3) 



with E(k) the energy spectrum such that the total energy 
is E = J E(k)dk. 

We examine two different flows, generated by differ- 
ent volume forces f that are either non-helical or fully 
helical (we consider a forcing function fully helical when 
(f • V x f ) is maximal and non-helical when (f ■ V x f) 
is zero, with the brackets denoting spatial average): 
the Taylor-Green (TG) vortex [23], and the Arn'old- 
Beltrami-Childress (ABC) flow [23|. The former is non- 
helical, and the resulting flow has no net helicity, al- 
though spatially localized regions with positive or neg- 
ative helicity develop. The latter is fully helical, and the 
resulting flow therefore has helicity (where the flow he- 
licity is defined as H = (u • V x u)). 

When using the TG vortex as a forcing function, we 
prescribe f as 

Ptg — fo [sra(kfx)cos(kfy)cos(kfz)x 

— cos(kfx) sm(kfy) cos(kfz)y] , (4) 

while the ABC forcing is given by 

Cabc = fo {[Bcos(k f y) + Csm(kfz)]x+ 
[Asin(kfx) + Ccos(fc/z)] y + 
[Acos(kfx) + Bsm(kfy)] z} . (5) 

Here fo is the forcing amplitude, which was set to have in 
the turbulent steady state all runs with r.m.s. velocities 
near unity. For ABC forcing, the constants were chosen 
to be A = 0.9, B = 1, and C = 1.1. Table [fl shows 
the parameters for all the runs. More details about the 
runs, and a detailed analysis of energy spectra, fluxes, 
and energy transfer, can be found in [a. |24| . 



III. THE SO(3) DECOMPOSITION 

The isotropic component of the longitudinal velocity 
structure functions for each flow is extracted using the 
SO (3) decomposition following the method described in 
[Tij . Some modifications were made considering the size 
of our simulations and to obtain a parallel version com- 
patible with the way data is distributed among comput- 
ing nodes. In this section we briefly introduce the method 
described in [Tij and give details of our implementation. 
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To do the SO (3) decomposition and recover the 
isotropic component, the longitudinal velocity structure 
function of order p is decomposed in terms of the spher- 
ical harmonics Yj m , obtaining 

S p (l) = ({[v(r)-v(r|l)].l} P ) 

= j2 s rm m (i), (6) 

where the brackets denote spatial average over the vari- 
able r, and homogeneity is assumed. The coefficients 
S° p m are obtained by projecting S p il) into the spherical 
harmonics, or equivalcntly, by projecting the p-th power 
of the longitudinal increments 

<Sw(r,l) = [v(r)-v(r + !)]•? (7) 

into the spherical harmonics and averaging over r. We 
are interested in the isotropic sector, given by the Sp°(l) 
functions, 

= ^ (J J £ SvP(r, l)P S m(<p)d<pd6^ , (8) 

where the angles ip and 9 are associated to the vector 1. 
The discrete version of this equation, and therefore the 
expression used in the numerical code to carry out the 
decomposition, is 

, N d N 3 

= ^ (9) 

j=l i=l 

where N is the linear resolution and Nd is the number of 
directions used to compute the average over the sphere. 

In Ref. [HI it was shown that 146 different directions lj 
covering in an approximately uniform way the sphere can 
be generated on a regular grid in such a way that all inte- 
ger multiples of lj lie on a grid point. This avoids the need 
to use three dimensional interpolations to compute the 
longitudinal increments Svfa, lj) when r + 1 does not lie 
on a grid point, significantly reducing the computational 
cost of the decomposition. The 146 directions are gen- 
erated by the vectors (1,0,0), (1,1,0), (1,1,1), (2,1,0), 
(2, 1, 1), (2, 2, 1), (3, 1, 0), (3, 1, 1) and those that are ob- 
tained by permuting their components in every possi- 
ble way (including multiplication of the vectors by — 1). 
With periodic boundary conditions, negative multiples of 
the vectors are not needed, since a spatial average over 
the entire box with the increment lj gives the same result 
as the average with —lj. This reduces Nd in Eq. © to 
73. 

The code we use to solve the Navier-Stokes equations 
is parallelized using a two-dimensional domain decompo- 
sition [2(| [ll| . Each computing node stores a slice of the 
velocity field in real space of size N x N x N z (N z < N, 
with N z a function of the number of computing nodes). 
As a result, increments in Eq. @ in the x and y direc- 
tions can be computed locally in each node. However, in- 
crements in the z direction require communication which 
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FIG. 1: Third order structure functions S3 (I) and isotropic 
component S3 (I) as function of / for the T3 run. Black curves 
correspond to the 73 different directions while the gray curve 
is the S3 (I) average. The straight lines indicate ~ / 3 and ~ I 
scaling. 

is handled using the Message Passing Interface (MPI) li- 
brary. The sum in Eq. (|9|) is then computed as follows: 
for each increment lj, displacements of the velocity field 
in the x and y direction are computed. Communication 
is then performed to displace the velocity field in the z 
direction if needed, and Sv(ji, lj) is computed for all val- 
ues of rj. Finally, the sum over all is done. The process 
is repeated for integer multiples of lj by just displacing 
the already displaced velocity field on lj over and over 
again. In this way, all communications are done between 
nearest neighbors avoiding all-to-all communications. 



IV. THE 4/5 LAW AND THE ENERGY 
SPECTRUM 

The result of computing the third order structure func- 
tions for the 73 directions using the T2 dataset is shown 
in Fig. [TJ The average over all directions (computed us- 
ing a procedure similar to the one described in [l4|) is 
also shown. From the 4/5 law J2{|, we expect the third 
order structure function to scale as the increment I in 
the inertial range. This is indeed the case for the average 
(the isotropic component) in a wide range of scales, while 
the structure function in each direction may or may not 
follow this law. In the dissipative range, where the flow 
becomes regular, a scaling ~ I 3 is observed for all curves. 

At lower Reynolds number (and spatial resolution) 
scaling in the inertial range can still be identified al- 
though the width of the scaling range decreases. This 
is illustrated in Fig. [3] with the isotropic second or- 
der structure function S$°(l) for runs Tl, T2, and T3. 
The scaling in the inertial range of this function is as- 
sociated with the scaling in the same range of scales of 
the isotropic energy spectrum. According to Kohnogorov 



4 



1 .000 








0.100 








0.010 


/ 


T2 


T1 


0.001 


T3 







0.01 0.10 1.00 

FIG. 2: Isotropic component of the second order structure 
function for runs Tl, T2, and T3 with increasing Reynolds 
number. The slopes of 2/3 and 2 are shown as references. 

theory (hereafter, K41) (25|, the second order structure 
function scales as ~ fa with (2 = 2/3, which implies in 
turn an energy spectrum E(k) ~ fc~^ 2+1 - ) = fc -5 / 3 . The 
slope of these structure functions in the inertial range 
(the range where S^°(l) ~ I) is slightly larger than 2/3, 
an effect associated to intermittency as discussed in more 
detail in the next section. As an example, a best fit to the 
power law in the second order structure function of run 
T3 gives C2 = 0.702 ± 0.004. Remarkably, at lower res- 
olutions we also recover deviations from the K41 predic- 
tion for £2 and the energy spectrum, with run T2 giving 
C 2 = 0.69 ± 0.01 and run Tl giving (2 = 0.68 ± 0.04. The 
second order isotropic structure functions for the ABC 
runs follow in the inertial range similar power laws, with 
slopes 0.66 ± 0.04, 0.695 ± 0.006 and 0.703 ± 0.003 re- 
spectively for the resolutions of 64 3 , 256 3 y 1024 3 grid 
points. The six values are within error bars. At large 
Reynolds numbers, S^CO ls c l° se to ~ l~ ' 7 and the en- 
ergy spectrum close to E(k) ~ k~ 1J for both helical and 
non-helical cases, which is slightly steeper than ~ fc~ 5 / 3 . 
This result is consistent with the result obtained from 
the largest simulation of isotropic and homogeneous tur- 
bulence done up to the moment using 4096 3 grid points 
fj, 26| , and with other simulations at large spatial reso- 
lution and Reynolds number (see e.g., Q). 

V. INTERMITTENCY 
A. High order moments and anomalous scaling 

Overall, the functions S P (V) for all values of p stud- 
ied display, after averaging over all directions, a dissipa- 
tive range that goes as ~ l p , an inertial range follow- 
ing some power law, and a range at large scales that 
depends on the forcing. Each individual direction be- 
haves as the average, although the structure functions 
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FIG. 3: Above: Sp°(l) structure functions as function of the 
increment I for p from 1 to 8 for the T3 run. S3 (I) is in- 
dicated by the thick curve. Only a range of scales near the 
inertial range is shown. Below: same for the isotropic struc- 
ture functions for run A3. 

for each direction show larger fluctuations and disper- 
sion, specially at large scales. This can be understood 
in terms of anisotropics associated to the forcing which 
prevail at large scales. The differences between partic- 
ular directions and the isotropic component decrease as 
the Reynolds number increases. 

According to K41 theory, the longitudinal structure 
functions of order p should scale in the inertial range 
as ~ fa with ( p = p/3, where ( p are the scaling ex- 
ponents. This scaling corresponds to a scale invariant 
(non-intermittent) flow. However, turbulence comes in 
gusts and regions with strong gradients do not fill the 
entire space. Strong events are localized in space and 
time, and the probability of finding strong gradients at 
small scales is larger than what can be expected from 
a Gaussian distribution. Intermittency leads to anoma- 
lous scaling (deviations from the ( p = p/3 relation), as 
the higher the order the more important are the contri- 
butions from these strong events. For the energy spec- 
trum, the scarcity of the small scales is responsible for 
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the steeper than k 5 / 3 energy spectrum discussed in the 
previous section. 

The dependence of the scaling exponent with the order 
is illustrated in Fig. [31 This figure shows the isotropic 
structure functions from p = 1 to 8 for runs T3 and A3. 
Only a range of scales near the inertial range is shown. 
In the inertial range, the longitudinal structure functions 
show scaling laws of the form Sp°(l) ~ as expected. 
As p increases, the scaling exponent £ p increases mono- 
tonically. To check convergence of the structure functions 
and exponents, we computed the accumulated moments 
0, Hj]> and verified that there is convergence up to the 
eight-order for these two runs (see e.g., [2J for the accu- 
mulated moments in run T3 for the average in only two 
particular directions). For the runs with lower resolu- 
tion, we computed structure functions up to the order p 
according to the convergence of their accumulated mo- 
ments. 

Let's now consider the scaling exponents £ p in the iner- 
tial range for all runs (where we define the inertial range 
as the range of scales where the 4/5 law holds). We 
start comparing the dependence of these exponents with 
the Reynolds number for the non-helical runs. Figure Q] 
shows £ p from the isotropic structure functions for the 
three runs with TG forcing. The f p = p/3 linear re- 
lation is indicated by the straight line. As previously 
mentioned, deviations of the exponents from the straight 
line are an indication of intermittency. Remarkably, the 
scaling exponents for the runs at lower Reynolds number 
are closer to the straight line than the exponents of the 
run at the largest Reynolds number, indicating intermit- 
tency still increases with the scale separation in the range 
of Reynolds numbers studied, and that more spatial res- 
olution is required to reach convergence with Reynolds 
number for high order moments. Figure 3] also shows the 
same results for the ABC runs. For this flow, a similar 
dependence of the exponents with the Reynolds number 
can be identified. 

Figure [5] shows the dependence of the second and 
fourth order scaling exponents for all runs as a func- 
tion of the Reynolds number. Because of the different 
forcing functions used, even when comparing runs at the 
same resolutions the Reynolds numbers arc slightly differ- 
ent. The second order exponent (2 slowly grows with Re, 
and seems to saturate for the largest Reynolds numbers 
reached. Note however that all values are within error 
bars. For (4 a decrease with the Reynolds number is ob- 
served, and in this case not all values are within error 
bars. The increase of the intermittency with Reynolds 
number is more pronounced when the intermittency ex- 
ponent /1 = 2(^3 — is studied as a function of Re (see 
Fig. only the runs at larger Reynolds numbers are 
shown, as determination of p, requires the sixth order ex- 
ponent). For the runs at the larger spatial resolution 
we obtain fi = 0.28 ± 0.02 and n = 0.30 ± 0.01 respec- 
tively for T3 and A3, while for runs T2 and A2 we find 
fx = 0.23 ± 0.07 and fi = 0.22 ± 0.03 respectively. Our 
values for the highest Reynolds numbers agree with mea- 





/ A - 




/ A 




/ t : 


- 


/■k - 







2 4 6 

P 



2.5 




2.0 


L / A - 

; A 4 ; 


1 .5 




1.0 




0.5 




0.0 





2 4 6 8 

P 

FIG. 4: Above: scaling exponents ( p in the T3 (triangles), 
T2 (crosses) and Tl runs (diamonds). The K41 prediction is 
given as reference. Below: same exponents for the A3 (trian- 
gles), A2 (crosses) and Al runs (diamonds). 



surements made by Meneveau and Sreenivasan [27J where 
they obtained /1 = 0.26 ± 0.03 (see also [28j]). Based in 
the results in Fig. [6l fx may be even larger for larger 
values of the Reynolds number. 



B. Role of Helicity 

To study whether previously reported differences be- 
tween helical and non-helical isotropic and homogeneous 
turbulence are associated to the presence of helicity or 
to anisotropics related to the forcing functions, we now 
focus on the datasets with the largest spatial resolution. 
Figure [7] shows the scaling exponents as a function of p 
up to the eight order for runs T3 and A3. The K41 pre- 
diction is shown as a reference, together with different 
models of intermittency: the log- normal model [29j | . 

Cp = | + ^p-p 2 )-, (10) 
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FIG. 5: Second order (above) and fourth order (below) scaling 
exponents as a function of the Reynolds number for all runs 
(triangles are for ABC forcing and crosses for TG forcing). 
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FIG. 6: Intermittency exponent as a function of the Reynolds 
number for the four runs at larger Reynolds number (labels 
as in Fig. [SJ). 



where fx is the intermittency exponent previously defined, 
the She-Leveque model [3G |. 



Cp _ 9 + 



1 
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the mean-field approximation [3 11 ] 
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p 3(l + 0.05p) J 
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FIG. 7: Comparison between the £ p exponents from the T3 
(crosses) and A3 runs (triangles). Shown as a reference are the 
K41 (solid atraigth line), log- normal (dash), mean field (dot- 
ted), She-Leveque (dash-dotted), and Arimitsu and Arimitsu 
(solid) predictions, the mean field model and Arimitsu and 
Arimitsu model are indistinguishable. The inset shows a zoom 
for the highest order moments with error bars. 



where 



Cp/3 - 1 



2(J) 2 (l- 9 )Ann(2), 



(14) 



(13) 



and the quantities cvo, X and q are determined from the 
intermittency exponent /i following the expressions in [32], 

m. 

The She-Leveque model and the mean field approxima- 
tion have no free parameters, while the log-normal model 
and the model of Arimitsu and Arimitsu depend solely on 
the intermittency exponent. For the log-normal model, 
we use n — 0.28 which is compatible with the value found 
in the A3 and T3 runs. Arimitsu and Arimitsu state that 
for infinite Reynolds number fi = 0.22 should be used in 
their model, and in [33| give explicit values for ceo, X and 
q for this case. These are the values we use in Fig. [7j 
and when plotting the model with values corresponding 
to larger values of \x we observed an improvement in the 
model although measurable differences for the highest or- 
ders persist. The best fit to both helical and non-helical 
data is given by the log-normal model, although it is well 
known that for higher order moments the model will fail 
as its exponents do not increase monotonically with p. 
The data deviates from all the other models. 

Table [HI shows a comparison of the £ p exponents of 
the isotropic sector for the helical and non-helical runs 
at the largest spatial resolution. Except for the highest 
order computed, differences are within error bars, and 
we can thus conclude that previously reported differences 
measured without the SO (3) decomposition were associ- 
ated to contributions from the anisotropic sector. The 
small discrepancy observed for p = 8 may be related to 
Reynolds number dependence, as runs A3 has a slightly 
larger Reynolds number than run T3, and as it was noted 
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TABLE II: Comparison of scaling exponents in helical and 
non-helical flows at the largest Reynolds number studied, p 
is the order, £a3 are the isotropic scaling exponents of the 
helical run, (t3 are the exponents of the non-helical run, and 
C-43 — Ct3 is their difference. 



CA3 



Ca3 — £? 



0.372 ± 0.002 
0.702 ± 0.004 
0.998 ± 0.005 
1.263 ± 0.007 
1.500 ± 0.009 
1.71 ± 0.01 
1.90 ±0.01 
2.07 ±0.01 



0.373 ± 0.001 
0.703 ± 0.003 
0.999 ± 0.004 
1.263 ± 0.005 
1.495 ± 0.006 
1.700 ± 0.007 
1.877 ± 0.008 
2.028 ± 0.009 



0.001 ± 0.003 
0.000 ± 0.007 
0.00 ±0.01 
-0.00 ±0.01 
-0.00 ±0.01 
-0.01 ±0.01 
-0.02 ±0.02 
-0.04 ±0.02 



before, the higher order exponents slowly decrease as the 
Reynolds number is increased. 



VI. CONCLUSIONS 

In this work we studied scaling exponents and intermit- 
tency in the isotropic sector of turbulent flows at differ- 
ent Reynolds numbers, generated by Taylor-Green forc- 
ing and by Arn'old-Beltrami-Childress forcing. The first 
is non-helical while the latter is maximally helical. The 
spatial resolutions were 64 3 , 256 3 and 1024 3 grid points, 
with the Reynolds number (based on the integral scale) 
ranging from « 40 to w 6200. 



A dependence of the intermittency corrections with the 
Reynolds number was found, with a larger intermittency 
exponent as the Reynolds number was increased. When 
the SO (3) decomposition was used, no dependence of the 
level of intermittency with the hclicity content in the 
flow was found, and small differences at the highest order 
studied may be due to Reynolds number dependence of 
the exponents. Finally, we confirm previous results in- 
dicating the SO (3) decomposition increases the range of 
scales where scaling is observed, and improves the deter- 
mination of scaling laws. 

When comparing with models of intermittency in the 
literature, it was found that the data is more intermit- 
tent than predictions from the She-Leveque model, the 
mean field approximation, and the model of Arimitsu 
and Arimitsu. In spite of the well known problems with 
the log-normal model, this model gives the best fit to 
our data. In this context, it is interesting to question 
the need to develop new models of intermittency with- 
out better data or a deeper understanding of the origin 
of intermittency in homogeneous and isotropic turbulent 
flows. 
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